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The quark condensate is calculated within the world-line effective-action formalism, 
by using for the Wilson loop an ansatz provided by the stochastic vacuum model. 
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CN . Starting with the relation between the quark and the gluon condensates in the heavy- 

\ quark limit, we diminish the current quark mass down to the value of the inverse 

■ vacuum correlation length, finding in this way a 64%-decrease in the absolute value 

of the quark condensate. In particular, we find that the conventional formula for the 
heavy-quark condensate cannot be applied to the c-quark, and that the corrections 
to this formula can reach 23% even in the case of the 6-quark. We also demonstrate 



that, for an exponential parametrization of the two-point correlation function of 
gluonic field strengths, the quark condensate does not depend on the non-confining 
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^ I non-perturbative interactions of the stochastic background Yang-Mills fields, 
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^ '. I. INTRODUCTION 
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As it is well known, the so-called chiral SUL(A^f) xSUR(A^f) symmetry of the classical action 
^ ■ of QCD with A^f massless quark flavors is spontaneously broken at the quantum level, with the 

i order parameter for this symmetry breaking being the quark condensate {ipip). Together with 

confinement, which is characterized by the gluon condensate {{gF'^^Y), chiral-symmetry breaking 
is one of the two most important non-perturbative phenomena in QCD. A natural question can 
be posed as whether these phenomena are interrelated or not. An affirmative answer to this 
question would imply proportionality between the quark and the gluon condensates, which was 
indeed found in Ref. The corresponding relation reads 

where Tg is the vacuum correlation length, at which the two-point gauge-invariant correlation 
function of gluonic field strengths exponentially falls off. Equation (fl]) stems from the integration 
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in the QCD partition function over soft gluonic fields in the leading, Gaussian, approximation. 
Within this approximation, the kernel of the four-quark interaction is defined by the two-point 
field-strength correlation function with the amplitude {{gF'i^^Y) and the correlation length Tg. 

Alternatively, if one first integrates in the QCD partition function over the quark fields, one 
arrives at a gauge-invariant effective action, where the gluonic degrees of freedom are represented 
in the form of Wilson loops and their correlation functions j2|. An advantageous feature of this 
approach is that, owing to the color- neutrality of Wilson loops, the calculation of the effective 
action becomes reduced to the calculation of the world-line integrals in an Abelian gauge theory. 
When the dynamical quarks which are integrated out are sufficiently heavy, namely their current 
mass M is larger than l/T^, the gluonic field inside the quark trajectory can be treated as nearly 
constant. In this heavy-quark limit, the one- loop effective action yields the following heavy-quark 
condensate of a given flavor jsl: 

^^^^^^^ = 

This expression coincides with the one known from the SVZ sum rules ^. 

With M decreasing downwards 1/Tg, variations of the gauge field inside the quark trajectory 
produce corrections to Eq. ([2]). The aim of the present paper is the calculation of such corrections. 
They will be obtained by using the approach of Ref. [5], which, for the case of a fermion moving 
in an arbitrary Abelian gauge field, yields a closed formula for the effective action with two field 
strengths. Furthermore, it is known that, in addition to the confining interactions of stochastic 
gluonic fields, there also exist non-confining non-perturbative interactions of those fields, albeit of 
a relatively small strength (cf. Ref. 6|). Below, we study the influence of such interactions on the 
heavy-quark condensate. For the case of the simplest, purely exponential, two-point correlation 
function of gluonic field strengths, we find the interesting phenomenon of a complete independence 
of the heavy-quark condensate from the non-confining non-perturbative interactions. 

The paper is organized as follows. In the next Section, we calculate the quark condensate 
by accounting in the effective action for the confining interactions of stochastic gluonic fields. In 
Section III, we generalize this result to the case where non-confining non-perturbative interactions 
of those fields are taken into account. Section IV provides a summary of the results obtained. 
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II. CORRECTIONS TO THE HEAVY-QUARK CONDENSATE 



Integrating over the quark fields in the QCD partition function, one arrives at the following 
one- loop effective action {2, 3, 5|: 



r°° ds 

{r[A;]) = 2NJ ^e- 
Jo s 
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Here, P and A stand, respectively, for the periodic and the antiperiodic boundary conditions, 
so that Jp = ^^^^^ = {s)=-ip (0)' ^ ^^'^ current quark mass. Since the quark 

condensate is always associated with a given flavor, we set iVf = 1. The corresponding expression 
for the quark condensate reads 
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where V is the Euclidean four-volume occupied by the system, and 
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is the average over gluonic fields. In the heavy-quark limit of M ^ ^/Tg-, the one-loop ap- 
proximation becomes exact, leading to Eq. ([2]) [cf. Ref. [Si] and the paragraph after Eq. (fT3|) 
below]. We notice that Eq. ([2]) uses the fact that the Yang-Mills field-strength tensor F^^ = 
d^A'j^ — dvA^^ — gf"'^'^A^^Al, which enters the quark spin term in the world-line action, can be 
recovered by means of the area-derivative operator acting on the Wilson loop j?!. By virtue 
of this fact, all the gauge-field dependence of the effective action becomes encoded in the Wilson 
loop. The latter is defined by the usual formula = {tiV exp (^ig JJ* rfr T^A^i^) ) , where 

is a generator of the SU(Ai"c)-group in the fundamental representation, and V denotes the path 
ordering. 

Since is completely determined by the geometric characteristics of the contour C = 

z^(r), the calculation of the quark condensate becomes an Abelian problem. In this section, 
we consider the confining part of {W{C)), deferring the study of the subleading non-perturbative 
non-confining part to the next section. Within the stochastic vacuum model ^, the corresponding 
area-dependent part of the Wilson loop reads 



{W{C)) = N, exp 
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In this formula, S is the minimal surface bounded by the contour C, and G = {{gF'j^^) ) is the gluon 



condensate. Furthermore, we choose the surface element (ia^j, in the form of an oriented, infinitely 
thin triangle built up of the position vector -2^(t) and the differential element dz^ = z^dr, namely 

)dT. One can then readily check that J da^y{z) J da^uiz') = (Jg* dr i^z^)^. 



da^u{z 



as it should be 
one has 



14|. Then, by virtue of an elementary Fourier transform J^e /^I^I+^p^ = .^-J^IL^^^^ 
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where E^,^ = T.^^{x;C) = ^ dr {z^Zi, — ZyZ^)5{x — z{t)) and f^ = f d^x. The exponential of 
interest thus reads 

where the antisymmetry of B^j^y has been used at the final step. One recognizes in this formula a 
Wilson loop corresponding to the Abelian gauge field 
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-z^B^y{z). 



The strength tensor of this field is F^^ = d^A^ - d^A^ = B^y + C^^, where C^u{z) = lzx{d^Bx,y - 
dvBxfji)- In particular, owing to just the Abelian Stokes' theorem, it is the strength tensor F^^ 
which automatically appears in the quark spin term of the effective action, being recovered by 
the operator -J— in Eq. ([3]). 



Accordingly, the one-loop effective action <^ takes the form 

ds e 



{V\Al]) = 2iV, 



F,y{x)mF^y{x) / , 
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with the corresponding Abelian covariant derivative Dn = — iA^ entering the formfactor 
F(^). In the spinor case at issue, this formfactor reads j5| F(^) = ^^^2^^ ~ \fiO^ where /(^) = 
Jp^ e"^^"""-*^ and ^ = sD'j^. In what follows, we find convenient to identically represent the 
formfactor F(^) in the form 



du 
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Following the method of Ref. [9], each of the two exponentials in the last expression can be 
represented as 
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At the final step of tlie transformation, by performing an elementary a-integration, we obtain for 
the formfactor F(^) a compact expression 

We now insert this expression into Eq. ([9]). Since the gauge field enters Eq. (fTOj) only via the 
exponential e^'"^'", we obtain for the 5-average in Eq. (|9]): 

/ F^,(a;)e^-^-F^,(a;)\ = V {F^MF,u{y))B = V [{B^,{Q)B^,{y))B + {B^MC^.u{y))B\. 

Jx / B 

(11) 

To obtain the last equality in Eq. f|TT]) we have used the fact that C^y{0) = [15]. The correlation 
functions that enter Eq. (fTTl) now read 



{B,UO)B,Ay))B = 12 j VBe'^^^''^-'"^'"^''"''B{Q)B{y) = ^f'^"^ (12) 

and 

{B,MCAy))B = ^ {yxd,) e-^l^l = ^^\y\ e-'^l^l (13) 

where we have taken into account that an antisymmetric tensor has 12 components. 

The large-M limit at issue corresponds to ,^ ^ 1 and ^\y\ ^ 1, so that f[^) — t- 1 + |, 
F(^) — )■ — |, and only Eq. (fT2|) contributes to Eq. (ITT]) , whereas Eq. (fT3l) does not. Recalling 
finally the definition of the quark condensate, Eq. (jlj), we recover Eq. ([2]). 

We will now apply the same method of calculation of the effective action to a derivation of the 
quark condensate for the smaller values of M, down to M = /i. Equations (fTU]) -f lT5]) yield 
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m^rrsy Jo Jy[y' 2Hl-uW)' 4 J' ^'^^ 

where "1" in the last bracket stems from Eq. (fT2l) . while ^— ^ j stems from Eq. (fT3|) . Accordingly, 
by using Eq. (jl]), we can obtain for the quark condensate the following expression: 
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The s-integration in this formula can be performed exactly. Denoting 

A 



and a 



we arrive at the following intermediate result: 
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where i^j^'s are the Macdonald functions. The 2;-integration here can still be performed analyti 
cally, yielding 



(^V')svz 

where (ijjip) svz is given by Eq. ([2]), and /(A) stands for the following integral: 
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For A ^ 1, the leading large-A terms 4 + (|)^ ■ "^fj^i yield /(A) — 1. For arbitrary A's, the 
remaining w-integration has been done numerically, with the result plotted in Fig. [H In particular, 
we obtain /(I) ^ 0.36, that is, a 64%-decrease in the value of the quark condensate when M ~ 
/i. In reality, only the values of J(A) corresponding to M = Mc, M;,, and Mt are of physical 
significance. We use the standard quark masses — 1.3 GeV, M/, ~ 4.2 GeV, and Mt — 173 GeV. 
The vacuum correlation length in full QCD with light flavors 
H ~ 580 MeV. This yields 



10] . Tg ~ 0.34 fm, corresponds to 



I{Mc/fJ,) ~ 0.60, /(MfcZ/i) ~ 0.84, /(Mt//i) ~ 0.996 in full QCD. 



(17) 



For the alternative case of quenched QCD, that is, SU(3) pure Yang-Mills theory, the vacuum 
correlation length is {ll| Tg ~ 0.22 fm, which corresponds to ~ 897 MeV. For this value of /i, 
we have 



I{Mjfi) ^ 0.47, I{Mb/fi) ~ 0.77, I{Mt/fi) - 0.993 in quenched QCD. 



(18) 



The sets of numbers (fT7j) and (fT8|) illustrate the degree of accuracy of Eq. ([2]) for various heavy 
flavors and various values of the vacuum correlation length Tg. Since the case of heavy quarks 
considered here constitutes an intermediate case between QCD with light quarks and quenched 
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Figure 1: The function /(A) in the range A G [1, 298], where 298 ~ ^ is the maximum value of A, which 
corresponds to Tg = 0.34 fm. 



QCD, the genuine value of I{Mf/^), for a given heavy flavor /, lies somewhere in between the two 
corresponding values of I{Mf/fi) listed in Eqs. (1171) and (ITS]) . In any case, we can conclude that 
Eq. ([2]) is inapplicable to the c-quark, since it can develop up to 53%-corrections [cf. I{MJa) 
from Eq. (fT8|) ]. We notice that a qualitatively similar conclusion has been drawn in Ref. 12|, 
where the leading correction to Eq. ([2]) has been evaluated through a non-perturbative gluon 
propagator in the Fock-Schwinger gauge. Finally, setting in Eq. ([2]) a certain heavy flavor /, and 
denoting ('?/''?/') svz,/ = — f^- ; '^^ can write, instead of Eq. f|T6|) . 



Mf 



(19) 



(V'^)svz,/ 

For an illustration, we plot in Fig. [2] the ratio (fT9|) for the case of / = 6 and Tg = 0.34 fm, up to 
M = Mb. In accordance with the intuitive expectations about the behavior of {ipip) with M, we 
observe a monotonic decrease of | {ipil^) \ with the increase of M. 
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Figure 2: Ratio (|19p for f = b and Tg = 0.34 fm, in the range of masses M G [fi, Mf,]. The corresponding 



maximum value of A is — ~ 7.24. 



III. ACCOUNTING FOR THE NON-CONFINING NON-PERTURBATIVE 

INTERACTIONS 



In addition to the confining interactions of stocfiastic background Yang-Mills fields, which lead 
to the Wilson loop in the form of Eq. , there also exist non-confining non-perturbative inter- 
actions of those fields. In this section, we demonstrate the interesting phenomenon of complete 
independence of the quark condensate from such interactions, provided they exhibit exponential 
correlations. 

To account for the non-confining non-perturbative interactions, one represents the full two- 
point correlation function of gluonic field strengths in the form jo, 8, 13 1 



Q §ab 1 — K 

Here, n G [0, 1] is some parameter, which defines the relative strength of the confining and the 
non-confining non-perturbative interactions. The lattice simulations in the SU(3) Yang-Mills 
theory yield the value of k = 0.83 ±0.03 (cf. Ref. [6i]), which means that the relative contribution 
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of the non-confining non-perturbative interactions amounts to only 17%. Expressing the Wilson 
loop via the correlation function {g'^F^^^{Qi)F\p{x)^ through the non-Abelian Stokes' theorem and 
the cumulant expansion [8], and using the above parametrization of (^g"^ F^^^^O) F^^^x)') , one obtains 
the following generalization of Eq. : 

(H/(C)>=iV.exp{-g^x 

K J da^^{x) j da^y{x') + ^—^^ j) dXfj, j) dx'^il + iJi\x - x'\) e"'"'^""''' |. (20) 
The non-confining non-perturbative interactions produce in the Wilson loop a term with the dou- 

Qoo 
I) §c '^^/^ §c '^•^M / ^'^ e~^"^ . 

The corresponding expression in Eq. (120|) resulted from the r-integration in this formula. 

Much as for the surface- dependent part of the Wilson loop, for the contour-dependent part 

we can also use some elementary Fourier transform, namely J^(l + /i|a;|)e~'^'^'"'"*^^' = ; to 

represent it as 

where = j^{x\ C) = dx^{T)6{x — x{t)). Further introducing a notation for the mean value 
(. . . )^ = / Vhf,e~^^^^^^^ ^^'"^'^'^^'^^^^'^^''^i- ■ ■ ), we notice that the full Wilson loop ([20]) can be 
written as a product of two averages: 

{W{C)) = (ehLB^.^^.^j . (^^i!^h^i^\j _ 

This equation generalizes Eq. (JTj) to the case where the non-confining non-perturbative interac- 
tions are also taken into account. Accordingly, the auxiliary Abelian gauge field ([S]) becomes 
now Ay{z) = ^z^Bfj_i,{z) + h^{z). Its strength tensor reads F^^, = Bfj_^ + Cfj,^ + H^i,, where 
Hiiu = d^hy — dyh^j_. Furthermore, Eq. (fTTl) also gets modified as 

! F^,(x)e^-^-F^,(x)\ = V ■ [{B^MB,u{y))B + {B^,{Q)C^M)b + (//^.(0)iJ^,(y)),] , 

/ B,h 

where we have denoted ((•■■) B)h as just (■ ■ ■ ) B,h- The appearing additional correlation function 
{H^y{{])H^j_v{y))h can be readily calculated by means of the formula 



10 



The result reads 



{H,MH,Ay))k = ^-^^^ ■ (/ibi -4) -e-^l^l. 



Using now Eqs. (fT2|) and (fT3|) . with G replaced by nG, we observe a remarkable mutual cancellation 
among all the Ac-dependent contributions. Namely, we obtain 



1 
V 



B,h 



Thus, Eq. (HM . with {■■■)b replaced by stays unchanged, and so does the resulting 

quark condensate. 

The question whether the obtained cancellation among the K-dependent contributions is spe- 
cific for the above-considered exponential ansatz for the correlation function (5f^F^j,(0)F^^(a;)), 
or it holds equally well for other ansatze (such as e.g. the Gaussian one), requires a separate 
study, which lies beyond the scope of the present paper. We only notice that, even in the absence 
of such a cancellation, the contribution of non-confining non-perturbative interactions is always 
suppressed, in comparison with the contribution of confining interactions, by a relative factor of 
^ ~ 0.2. 



IV. SUMMARY 



The aim of the present paper was to find a relation between the quark and the gluon con- 
densates, which would yield, for various heavy fiavors, corrections to the known Eq. ([2]). The 
corrections thus obtained, given by Eqs. (fT7j) and (fTSl) . show that Eq. ([2]) applies with a good 
accuracy only to the t-quark. Rather, for the 6-quark, the corrections are ~ 20%, while for the 
c-quark they can be as large as ~ 50%, thereby making Eq. inapplicable to the c- and the 
s-quarks. Also, as one can see from Fig. 1, when the continuously varied current quark mass 
M reaches the value of the inverse vacuum correlation length /i, the absolute value of the quark 
condensate decreases by 64% compared to the value provided by Eq. ([2]). 

We have used in our calculations the most general ansatz for the Wilson loop, which is pro- 
vided by the stochastic vacuum model and accounts for the confining and non-perturbative non- 
confining interactions of the stochastic gluonic fields. The corresponding two-point surface-surface 
and contour-contour self-interactions of the Wilson loop can be represented as being mediated 
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by an auxiliary Abelian gauge field with the Gaussian action. In particular, for the most simple, 
exponential, parametrization of the two-point correlation function of gluonic field strengths, we 
have found an interesting phenomenon of a complete independence of the heavy-quark condensate 
from the non-confining non-perturbative interactions of the stochastic gluonic fields. 

In conclusion, we have started our analysis from the heavy-quark limit, where chiral sym- 
metry is explicitly broken by a large current quark mass. The advantage of working in this 
limit is that one avoids possible uncertainties related to the particular form of the field-strength 
correlation function. Indeed, owing to the constancy of the gauge field inside the heavy-quark 
trajectory, Eq. ([2]) in the t-quark case turns out to be almost exact. We emphasize that even 
in the heavy-quark limit we still have a relation connecting the quark condensate {ipip) with the 
gluon condensate {{gF^j^Y'). We have not proceeded to the current quark masses smaller than 
(cf. Fig. 1), which is the case of s-, c?-, and -u-quarks. The reason is that, for such light quarks, the 
effect of spontaneous breaking of chiral symmetry starts to play an important role, resulting in 
the appearance of a significant self-energy contribution to the dynamical constituent quark mass. 
Thus, since such a self-energy contribution cannot be consistently calculated within the adopted 
world-line formalism, we have to restrict our analysis to the case of heavy quarks, for which this 
contribution can be safely disregarded compared to the current quark mass. However, even in 
the heavy-quark case provided by the h- and c-quarks, we have found substantial corrections to 
Eq. ([2]). The way in which Eq. ([2]) along with these corrections goes over into Eq. ([T]) for light 
quarks can be the subject of a separate study. 
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